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In this paper we discuss the properties of the reduced density matrix of quantum many body 
systems with permutational symmetry and present basic quantification of the entanglement in terms 
of the von Neumann (VNE), Renyi and Tsallis entropies. In particular, we show, on the specific 
example of the spin 1/2 Heisenberg model, how the RDM acquires a block diagonal form with 
respect to the quantum number k fixing the polarization in the subsystem conservation of S z and 
with respect to the irreducible representations of the S n group. Analytical expression for the RDM 
elements and for the RDM spectrum are derived for states of arbitrary permutational symmetry 
and for arbitrary polarizations. The temperature dependence and scaling of the VNE across a 
finite temperature phase transition is discussed and the RDM moments and the Renyi and Tsallis 
entropies calculated both for symmetric ground states of the Heisenberg chain and for maximally 
mixed states. 
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I. INTRODUCTION 

Entanglement properties of interacting quantum 
many-body systems [I| lies at the heart of many quantum 
information processes such as measurement based quan- 
tum computation, teleportation, security of quantum key 
distribution protocols, super-dense coding, etc. 0. Be- 
ing a principal resource for quantum information, one is 
interested to know how much entanglement is present in a 
system and how much of it can be used or created. Entan- 
glement also provides benchmarks for success of quantum 
experiments. Entanglement properties are presently in- 
vestigated for several spin chains Q, Q, PjB, 



for strong! 



and pairing models 



yc< 

tei 



correlated fcrmions [6|, 
for itinerant bosons 



13] 
etc. 

The calculation of the entanglement involves the 
knowledge of the reduced density matrix (RDM) char- 
acterizing quantum systems in contact with the envi- 
ronment such as a thermal bath or a larger system of 
which the original system is a subsystem. In particu- 
lar, the spectrum of the RDM, which by definition is 
real and nonncgative with all eigenvalues summing up 
to one, provides an intrinsic characterization of a sub- 
system. The relative importance of a subsystem state, 
indeed, is directly related to the weight that the corre- 
sponding eigenvalue has in the RDM spectrum. Thus, 
for instance, the fact that the eigenvalues Xi of the RDM 
for a one dimensional quantum interacting subsystems 
decay exponentially with i implies that the properties of 
the subsystems arc determined by only a few states. This 
property is crucial for the success of the density-matrix 
renormalization group (DMRG) method [l(| in one di- 
mension. In two dimensions this property is lost [l7| and 
the DMRG method fails. 

For a subsystem consisting of n sites (or n q-bits) the 



RDM is of rank 2™ so that for large n the calculation of 
the spectrum becomes a problem of exponential difficulty 
While the spectrum of the full RDM for subsystems with 
a small number of sites (e.g. n < 6) has been calculated 
[lH, the full RDM for arbitrary n is, to our knowledge, 
exactly known only for the very special case of non inter- 
acting quantum systems such as free fermions (see e.g. 
[l9j ) or free bosons. 

The aim of the present paper is to analytically calcu- 
late the elements of the RDM of permutational invariant 
quantum systems of arbitrary size L, for arbitrary per- 
mutational symmetry of the state of the system (labelled 
by an integer number < r < L/2) and arbitrary sizes n 
(number of q-bits) of the subsystem. We remark that the 
invariance under the permutational group physically im- 
plies that the interactions among sites have infinite range. 
As an example of such system we consider the Heisen- 
berg model of spin 1/2 on a full graph consisting of L 
sites, with fixed value of spin polarization 5* 2 = L/2 — N. 
For this system we calculate the RDM and the entangle- 
ment von Neumann entropy (VNE) for a subsystem of 
arbitrary n > 1 sites for arbitrary L, N. The tempera- 
ture dependence and the scaling properties of the VNE 
across a finite temperature phase transition occurring in 
the system are also discussed, and the RDM moments 
and the Renyi and Tsallis entropies calculated both for 
symmetric ground states of the Heisenberg chain and for 
maximally mixed states. 

The plan of the paper is the following. In Section U we 
discuss model equations and provide basic definitions. In 
Section |TT] we consider the main properties of RDM el- 
ements and show how the symmetry properties of the 
system allow to decompose the RDM into a block diag- 
onal form. In Sec. IIVI we present an exact analytical 
expression of the RDM matrix elements for arbitrary pa- 
rameters values whose rigorous proof is provided in the 
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thermodynamic limit in the appendix [B] and for the case 
of fully symmetric states in the appendix fX] In Sec. [V] 
we provide an analytical characterization of the RDM 
spectrum and discuss scaling properties and tempera- 
ture dependence of the von Neumann entropy. In Sec. 
VII moments of the reduced density matrix are discussed 
and several quantities of interest like mutual information, 
Renyi and Tsallis entropies, are calculated. In the last 
section the main results of the paper arc briefly summa- 
rized. 



II. MODEL EQUATION AND BASIC 
DEFINITIONS 



We consider a permutational invariant system of L 
spins 1/2 on a complete graph with fixed total spin polar- 
ization S z = L/2 — N and described by the Hamiltonian 



2L\ 2 \ 2 



1 



hS, 



(1) 



Here S = (S x 



Sy, S Z 



S a = iEf=i^?> with of Pauli 



matrices acting on the factorized Y[ ®C2 space. This 

l 

Hamiltonian is invariant under the action of the sym- 
metric group Sl and conserves S z , [H,S Z ] =0. A 
complete set of eigenstates of H are states \^ l.n.t) asso- 
ciated to filled Young tableaux (YT) of type {L — r, r}(N) 
(see [2l[ for details), where the subscript N denotes the 
number of quanta present in the tableau and the symbol 
{L — r, r} refers to a tableau of only two rows, with L — r 
boxes (sites) in the first row and r in the second row. For 
the Hamiltonian in ([1} we have: 



H^L^,)-) — EL,N,r\^L,N,r), 

1 fJr 



El,n,t — - 

S z \^L,N,r) 



L 



(L - r + 1) + h(L - 2N) 



(2) 



L 



N)\y L ,N.r) 



where TV = 0,1,..., L determines possible values 
of the spin polarization and r takes values r = 
0, 1, max(iV, L — N). Notice that, due to the sym- 
metry and antisymmetry of a YT with respect to rows 
and columns, respectively, the state \^L,N,r) can exist 
only if N > r (for the explicit form of the state \^L,N,r) 
see Eq. (|B3|) below). The degeneracies of the eigenvalues 
El,n,t are given by the dimension of the corresponding 
YTs: 



deg L . r = 



L 

r - 1 



(3) 



Consider a set of vectors \^ u ), u — 1, ... dcg L r , forming 
an orthonormal basis in the eigenspace of H with eigen- 
value El n r- We define the density matrix of the whole 



system as 



deg z 



E 



(4) 



It can be easily shown that UL,N,r possess the following 
properties: 

i) The matrix <JL t N,r has eigenvalues Ai = A2 = ... = 
Adeg L r = (deg£ r ) -1 , with remaining 2 L — dcg £ r eigen- 
values all equal to zero. This follows from the fact that 
each vector Y& u ) is an eigenvector of <?l.n .r with eigen- 
value 



Sl. 



Since the spectrum of (JL,N,r is real and 
nonnegative with all eigenvalues summing up to 1, the 
remaining 2 L — deg L r eigenvalues must vanish. 



ii) Matrix ul,n,t satisfies: {(JL,N,rY 



dcg £ 



-&L.N.r 



This 



follows from the definition (j4j and the orthonormality 
condition (*„,|* tt ) = 6 UW . 

iii) Matrices &L,N,r commute with each other 
[ci,iV,r i &L,N',r'] = 0. This follows from orthogo- 
nality of eigenspaces of H for different eigenvalues. 

iv) Introducing the operator Py, permuting subspaces i 

L 

and j of the Hilbert space Y[ on which the matrix 

1 

ol.n.t acts, we have that: [a, Py] = for any i,j. 
This last property can be proved by considering 



PijO-L,N,rPij = -, P ijl*t.)(*«|-f5 



V 

utJ BZ,,r u=1 



(5) 



The vectors |*^} = P 
being (ttJJO = ^ w \P^P i3 \^ u 



ij\^ u ) form an orthonormal basis, 



I*, 



= 6. 



WW : 



and (Pij) 2 = I- Now, the sum 



because = P^ , 

Y,t=i' r |*u)(*u| = ^dc gL ,. is a unity operator in a fac- 
tor space of dimension deg Lr , and therefore it does not 
depend on the choice of the basis. Note that vector \^f' u ) 
belongs to the same factor space as \^ u ), because per- 
mutation only results in different enumeration. Con- 
sequently, 



Pij ° P ij 



p, or [a, P^] = 0. 



(6) 



The latter property implies that in Eq. (HJ) the sum over 
the orthogonalizcd set of basis vector in (0} can be re- 
placed by the symmetrization of the density matrix di- 
rectly, namely a L ,N,r = ±\ J2p |*12...l)(*12...z,|) where 
the sum is over all L\ permutations of indexes 1,2, ...L, 
and \^i2...l) is some unit eigenvector of H with eigen- 
value pL,Af,r- In particular, it is convenient to choose 

|*12.. .L) = \#L,N,r), 



<7L,N,r = J\ ^ \^L,N,r)(^L,N,r\ 
P 



(7) 
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It is evident that such a sum is invariant with respect 
to permutations and that UL,N,r is properly normalized: 
Tra L . N . r = 1. 

The Reduced Density Matrix (RDM) of a subsystem of 
n sites is defined by tracing out L — n degrees of freedom 
from the density matrix of the whole system: 

P(n) = Tr L - n (JL,N,r- (8) 

Due to the properties © and |(5J), p/ n ) does not depend 
on the particular choice of the n sites, and satisfies the 
property §6§ in its subspace (we omit the explicit depen- 
dence of pr n \ on L, N, r for brevity of notations). 

III. RDM PROPERTIES AND BLOCK 
DIAGONAL FORM 

The RDM can be calculated in the natural basis by us- 
ing its definition in terms of observables: (/) = Tr(pt n )f) 
where / is a physical operator acting on the Hilbcrt space 
of the 2 n x 2" subsystem. The knowledge of the full set 
of observables determines the RDM uniquely. Indeed, if 
we introduce the natural basis in the Hilbert space of the 

n 

subsystem, Y[ ®C 2 , the elements of the RDM in this 

k=l 

basis are 

4tt = <^tt> = Tr (^)^tt)' ( 9 ) 

with = II ®z) k k and ej- a 2 x 2 matrix with 

1 2 "' " k=l k 

elements (e}) = SikSji. The matrix e^*?""^ has only 
one nonzero element, equal to 1, at the crossing of the 
row 2"- 1 ii + 2 ,l ~ 2 i 2 + ... + i„ + l and the column 2 n " 1 j 1 + 
2 n ~ 2 j2 + ■•■ + jn + 1 (all indices i,j take binary values 
i k = 0,1 and j k = 0,1). To determine all the RDM 
elements one must find a complete set of observables and 
compute the averages (e£*?' "^)- Note that a generic 
property of the RDM elements, which follows directly 
from (JB1) , is that any permutation between pairs of indices 
(im,jm) and (ik,jk) does not change its value, e.g. 

Jli2-in =p i ? il-in = „V>- i ? = ... = flM"- 1 "* . (10) 

Another property of the RDM follows from the S z invari- 
ance 

l>/;r', = 0, if H + *2 + + In ^31+ 31 + ... + jn- 

(11) 

Thus, for instance, the RDM for n — 2 has only 6 nonzero 
(4 different) elements , Poo;Poi = p%, p\l = p%\ and p\\, 
subject to normalization Trpoz) — Poo+ Pio + Pol + Pii = 
1. It is convenient to introduce the operators 

e = ^ Q J =a , e 1= ^ Q j = <r , 

e «= [o o) = p > e ^{o i)= h - 



If we represent a site spin up with the vector (J) and a site 
spin down with the vector ( ( j ) ) then p& and /ifc are spin up 
and spin down number operators on site k, while cj~ , a + 
represent spin lowering and rising operators, respectively. 
Thus, for instance, the observable (pij52^3^4-..^n) = 
Poo li "i gi ves th e probability to find spins down at sites 
3,4, ...n, and spins up at sites 1,2, while the observable 

( tT i" (T 2" cr 3~ (T 4~^5---^-n) = Piiooi!!!i gi yes the probability to 
find spins down at sites 5,6, ...n, spin lowering at sites 
3,4 and spin rising at sites 1,2. Note that the latter 
operator conserves the total spin polarization since the 
number of lowering and rising operators is the same. Also 
note that the correlation functions with a non conserved 
polarization vanish, e.g. 

(a+a+^h 4 h 5 ...h n ) = pTilllzl = 0, (12) 

in accordance with (fTl"]) . 

One can take advantage of the S z invariance (e.g. 
Eq. ITTj) to block diagonalize the RDM into independent 
blocks Bk of fixed polarization k = i\ + i% + ... + i n = 
ji + 32 + • • • + jn (here k = 0, . . . , n gives the number of spin 
up present in the subsystem) . In Fig. Q] the blocks Bk 
appearing in the RDM p(K\ have been shown for the case 
n = 5, L = 18, N = 8, r = 6. We remark that the n + 1 
diagonal blocks correspond to the values s z — (n— 2fc)/2, 
k = 0, 1, n the subsystem polarization can assume, be- 
ing the block decomposition a direct consequence of the 
S z symmetry. The dimension of the block Bk coincides 
therefore with the number of possible configurations that 
k spin up can assume on n sites, e.g. dimBfc = (V). One 
can check that the sum of the dimensions of all blocks 
gives the full RDM dimension, i.e. J2 k dimBk = 2". No- 
tice that the block diagonal form in the natural basis is 
achieved only after a number of permutations of rows and 
columns of the RDM have been performed. We also re- 
mark that the fact that the middle block B% in Fig. [T]has 
all vanishing anti-diagonal elements is purely accidental 
(see also remark at the end of Sec lIT|) . 

Blocks B k consist of elements eV ® e* 2 ® el 3 <g> . . . ® e\ n of 

J i. ,72 J 3 J ti 

n n 

the original matrix, with "^ip = "^2 j P ~ k and i p = 0,1, 

l l 

jp = 0, 1. In its turn, all elements <g> e 1 ^ (8) <g> .. . ® ej" 
of the block B k can be further block diagonalized (see 
below). In the natural basis, this diagonalization is done 
according to the number of pairs of type (e^CgieJ) present 
in the elements. In the following we denote by Gz the 
part of the block associated to elements with Z pairs 
(e?<8>eo) in it. The sub-block Go of the block B k is formed 
by the elements containing k terms e\ and (n — k) terms 
e[] in the product, i.e. e\ ® ...e\ ® e^® ... ® and all 
permutations. All such elements lie on the diagonal, and 
vice versa, each diagonal element of B k belongs to Go- 
Consequently, the sub-block Go consists of (?) elements. 
The number of elements, degGi(fc), in the sub-block G\ 
is equal to the number of elements of the type <8> ej <S> 
e-J(g)e- 2 (g)...(8)e^^, such that l+0+i 1 +i 2 +...+i n -2 = k. 
Using elementary combinatorics we obtain: 
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1 10 20 32 1 10 20 32 




1 10 20 32 1 10 20 32 



Figure 1: Left panel. RDM prg) in the natural basis. Param- 
eters are L = 18, N = 8, r = 6. Black boxes denote non zero 
elements. Right panel. The same matrix of the left panel af- 
ter the following chain of permutations of rows and columns 
have been applied to it: 7?8,i4-Ri4, 18-^8,21^3,29^2, 28^6,16^4, 9. 
(Ri,j denote the operator which exchange first columns i and 
j and then rows i and j). Black boxes denote non zero el- 
ements. Block diagonal structure associated with values of 
k = 0, 1, 5 is evident. The single element present in blocks 
k = 0, 5, are 1/34 and 1/153, respectively. Elements values 
inside other k-blocks are given in Fig. 

~=(DQ(rD- ™ 

Analogous calculations for arbitrary sub-block Gz yields 

From the restriction JZi *p = 52™ Jp = k we de- 
duce that the block Bk contains non-empty parts 
Go, G\, G f m in(fc,n-fe), leading to the following decom- 
position: 

min(fc,n— k) 

B k = [j G z . (15) 
z=o 

Indeed, the normalization condition following from (1151) . 
gives 

min(fe,n— fc) / \ 2 

£ degG z (*)= (") . (16) 

It is important to note that all elements of Gz are equal. 
This is a direct consequence of the property in Eq. (fTUf . 
A graphical representation of the Bk block for a partic- 
ular choice of L, N, r, is given in Fig(2] 

It is instructive to discuss the structure of blocks B k 
in terms of the matrices a in in ([7]) since this structure is 
directly connected with the block diagonalization of the 
RDM with respect the the irreps of S n . We have, indeed, 
that each block Bk can be decomposed in the form 

min(fe,n— k) 
s=Q 



1 2 3 4 5 12345678910 12345678910 




1 2 3 4 5 12345678910 12345678910 



Figure 2: Blocks for k = 1, 4 (left panel), k = 2 (central 
panel) and k = 3 (right panel) corresponding to the block 
diagonal form of p(5) in Fig[T] In the left panel white and 
black boxes denote elements 5/306 and 1/2448 for k = 1 and 
elements 2/51 and 1/1020 for case k — 4. In the center panel 
white boxes denote the element 1/1020 while black boxes, ex- 
cept for the ones along the diagonal which are all equal to 
2/51, denote the element —1/1190. In the right panel white 
boxes denote the element 1/1360 while black boxes denote the 
element —3/4760, except the ones along the diagonal which 
are equal to 1/34. The matrices cr n ,k,s involved in the decom- 
position in (|19|) have the same shapes as blocks Bk in the left 
panel for 05 n = (7541, with white and black boxes denoting el- 
ements 1/5 and —1/20, respectively. The shape of matrix 17521 
(resp. CT522) is the same as in the central panel, with white 
boxes denoting 1/60 (resp. —1/30 ) and black ones —1/15 
(resp. 1/30 ), except the ones along the diagonal which are 
1/10, and shape of (7531 (resp. 0532) is the shape is as in right 
panel with white box denoting 1/60 (resp. —1/30) and black 
ones —1/15 (resp. 1/30 ) except for the ones along the diag- 
onal which are given by 1/10. The 5 x 5matrices (7510 = (7540, 
and 10 x 10 matrices (7520 = 0-530 have all elements equal to 
1/5 and 1/10, respectively, while (7500 = (7550 = 1. 

where cr n ,k,s are associated to filled YTs of type {n — 
s, s}(fc) and a^'^'J are coefficients related to the corre- 
sponding eigenvalues X^k's °f * ne RDM" by 

L,N,r \ L,N,r j /-, \ 

a n,k,s = K,k,s de Sn,s ■ (18) 

Notice that the matrices <7 n .k.s in the natural basis 
have dimension 2" and coincide with the ones given in 
(U]) . In the proper basis (e.g. that of the irrep of S n ) 
they have dimension deg„ s x deg„ s and contribute to 
Bk with a sub-block of dimension deg n s corresponding 
to the filled YT of type {n — s, s}(fc). In performing this 
reduction one actually achieves the diagonalization of the 
block Bk, as it is evident from (fl~7|) (recall that <J n .k.s have 
eigenvalues 1/ deg„ s ). A first reduction of the matrices 
is achieved by accounting for the S z symmetry discussed 
before, this leading to matrices cr n ,fc,s of size (^) x (^). 
In Fig. [2] are also given the matrices o~ n ,k,s appearing 
in the decomposition of blocks Bk for the specific exam- 
ple considered in Fig. [T] and for k — 0, ...5, s = 0, 1, 2. 
One can check that these matrices satisfy all properties 
of matrix a given above and in particular, the number 
of their nonzero eigenvalues (all equal to 1/ deg„ s ) coin- 
cides with the dimension deg„ s of the YT to which they 
are associated. This implies that they can be further 
reduced from (£) x to deg n s x deg„ s size by elimi- 
nating the spurious (£) — deg„ s zero eigenvalues (these 
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eigenvalues arise because in the natural basis the dimen- 
sion of the representation is larger than the one of the 
S n irrcps). This is achieved by using the singular valued 
decomposition of the matrix a to write it in the form: 
a = UWV T , where W a diagonal matrix whose elements 
are the singular values and U and V are orthogonal ma- 
trices: U T U = V T V = 1, with superscript T denoting 
the transpose (this decomposition can be obtained very 
efficiently numerically (22j). 

The reduction to the sub-blocks of Bk in the proper S n 
representation is then achieved as: a = uwv T where u 
and v are rectangular matrices of dimension (?) x deg„ s 
obtained from U and V by omitting the columns corre- 
sponding to the zero eigenvalues and the matrix w is a 
deg ra s x deg„ s diagonal matrix with the nonzero eigen- 
values along the diagonal (in our case, since the nonzero 
eigenvalues of a are all equal to 1, w reduces to an unit 
matrix). The matrix w then provides the representation 
of u n ,k,s in the proper S n space leading to the full diag- 
onalization of the block Bf.. 

Thus, for example, the block diagonal form of the RDM 
in the right panels of Fig. Q] (see also Fig|2]) is expressed 
in terms of matrices a as 

B\ = A510O510 © 4A511 = A540CT540 © 4 A541CT541 = -B4, 
B 2 = A520 0-520 © 4 A521 0521 © 5 XA5220522, (19) 
B3 = A530 CT530 © 4 A531 0-531 © 5 x A532 CT532, 

(20) 



with blocks B = A 500 = 53, B 5 = A 550 = and 
eigenvalues A510 = 11/612, A 52 o = 76/1785, A530 = 
19/595, A540 = 11/255, of degeneracy 1, eigenvalues 
A511 = 13/816, A521 = 299/7140, A 53 i = 299/9520, 
A541 = 13/340, of degeneracy 4, and eigenvalues A522 = 
13/357, A532 = 13/476, of degeneracy 5 (having adopted 
the short notation X n k s = A^'^ r , the chosen parameters 
L = 18, N = 8, r = 6 are understood). 



IV. ANALYTICAL EXPRESSION OF RDM 
ELEMENTS 

The main analytical property of the RDM is summa- 
rized in the following statement: 

Elements gz of a sub-block Gz of a block Bk of the RDM 
for arbitrary L, N, r, n, are given by: 



gz 



IL-n-s 
\N-k) 



Z 



)( N z)( 1 



(-l)' r 



m=0 



N — r 
Z — m 



L-N- 
Z — m J \ rr 
(21) 



This expression has been derived by extrapolating exact 
results obtained for finite size calculations using symbolic 
programs and its correctness has been checked by com- 
paring with brute force numerical calculation of the RDM 
up to large sizes. Notice that Eq. (f2Tj) completely defines 
all elements of the RDM in the natural basis. In prac- 
tice, to find the element P, Q of the RDM (p( n )) P Q m the 



natural basis one must take the binary representation of 
numbers P — 1 and Q — 1 (which provide the sets of inte- 
gers {i p } and {j p }, respectively), find the corresponding 
number Z and use (|21[) . 

A proof of the statement for arbitrary L,N,n is given 
in Appendix A for the specific case r = corresponding 
to fully symmetric states. A proof of Eq. (|2Tj) which is 
valid in the thermodynamical limit L — > 00 is provided 
in Appendix B. In this respect, we remark that in the 
limit L — > 00 Eq. ((2"Tj) simplifies to 



9z =P 

where we denote with p 



i — k 



k„Z 



(i-pVv 



^ u 

L ' P 



and 



V 



(p-m)(i - p-m) 



p(l - p) 

For a proof of Eq (|22|) see Appendix [ 



(22) 



(23) 



V. SPECTRAL PROPERTIES OF RDM AND 
ENTANGLEMENT ENTROPY 

The existence of two representations for the block Bk 
of the RDM, one in terms of matrices Gz given in Eq. 
(|15| . the other involving matrices a n ,k,s and given in Eq. 
(TIT)) , have been shown in SccHIl These representations, 
together with the invariancc of Gz and o~ n ^k,s with re- 
spect to permutations, imply the existence of linear rela- 
tions of the form 



min(fc,n— k) 



G z {k) 



^2 Pz{k,s)a n 



(24) 



s=0 



where (3z{k, s) are constants and Gz denotes the matrix 
formed by all elements of Gz- Since o„ ^ s commute for 
different s (see the property (iii) of matrices a in Sec. 
we have that also Gz commute 



[G z (k),G z >(k)} = 0. 



(25) 



This also implies that all RDM eigenvalues \ n ' k 'J must 
be linear combinations of elements gz of matrices Gz- 
One can show, indeed, that the general expression of the 
RDM eigenvalues is 



, L,N'. 



min(k.n-k) 

E 

z=o 



Uz\n, k)g z , 



(26) 



with coefficients q^' (n, k) given by 



k-t 



.,0) 



i=0 



k — s\ (n — k + s 



Z-i 
(27) 

where Z,s = 0, 1, k. From Eq. (|2"T|) one can see that 

(n, k) are integer coefficients which, due to the prop- 
erty (f2"S")) . do not depend on the characteristics of the orig- 
inal state L, N and r. Thus the dependence of the RDM 



6 



eigenvalues on these parameters enters only through the 
elements gz PTjl . Moreover, one can shown that they 
satisfy the following relations 



(n,fc) = 1, 



(28) 
(29) 

(30) 



min(fc.n — k) 
min(fc,n— k) 

a { z\n,k)=0 for s > 0, 

min(fc,n— fc) f 7\ 

Y a^(n,k)( = 0,p = 0, (31) 

min(fc,n-fc) . . 

z=o 



<W) = (-1)\ 

(32) 



a proof of which can be found for special cases in [25|. 

From Eqs. ([2"T]l. ([22]). ([2"rj]l. (|2"7)) . the explicit analytical 
form of the complete spectrum of the RDM is obtained. 

The knowledge of the RDM spectrum allows to inves- 
tigate the bipartite entanglement, e.g. the entanglement 
of a subsystem of size n with respect to the rest of the 
system (see [l| for a review). This is done in terms of 
the entanglement entropy which for pure states at zero 
temperature coincides with the non Neumann entropy 

S( n ) = -tr(p {n) log 2 p (n) ) = - E Afc log 2 (33) 

where Afc the eigenvalues of the RDM p(„), obtained 
from the density matrix p of the whole system as p^ = 
t- r (L-n)P- For the infinite range ferromagnetic Heiscn- 
berg model at zero temperature the density matrix of 
the whole system is a projector on the symmetric ground 
state p = \^(L,N))(^(L,N)\ considered in [h| where 
it was shown that A fc = (fc)(jv~ fc)/(jv), where k = 
0, 1, ... min(n, N). In the limit of large n the VNE be- 
comes 



.5', 



1 



1 



n(L — n) 



- 2 lo S2( 27re M) + \ lo S2 ^ 



(34) 



One can show that a zero temperature (e.g. p = r/L = 0) 
[loj the spectrum of the reduced density matrix is de- 
scribed by a binomial distribution Afc 
converges to a Gaussian for large n 



pV" _fe (fc) which 



k \2 



A fc 



V27 



: CXp 



2a 2 



(35) 



where a 2 = np(l — p) 1. For finite temperature one 
introduces the thermal VNE for a block of size n as 

S(n)(fi) = -Hhn) W) l0 §2 P(n) 08)), (36) 
hn){P) = ^Y, d r e ~ PEr (P{n){r)), (37) 



0.25 



0.20 - 



0.15 



0.10 - 



0.05 - 



0.00 




Figure 3: Traces of the powers a = 2, 2.5, 3, 3.5, 4 (from top 
to bottom, respectively) of the RDM versus the block size n 
for symmetric states r = 0, as obtained from the analytical 
expression (|47p for p = 0.5 (continuous curves) and from exact 
expressions of RDM eigenvalues (|26[) . (|27[) (full dots) for L = 
2 * 10 3 . 



where p(„)(/3) is the thermal reduced density matrix, Z 
is the partition function, and (.) denotes the equal weight 
(thermic) average over all orthogonal degenerate states, 
corresponding to a given permutational symmetry. Note 
that (p( ra )( r )) commutes with any permutation opera- 
tor and does not depend on the choice of sites in the 
block but only on its size n. Also note that the matrices 
(p( n )( r )) commute for different r 



[(P(n)(r)),(p(n)(r'))] =0, 



(38) 



so that the diagonalization of p^((3) is reduced to the 
diagonalization of (p(n)( r )) f° r arbitrary r. From Eq. 
P7[) the computation of the temperature-dependent von 
Neumann entropy is easily made with the help of the 
general expression of the eigenvalues of the RDM in Eq.s 
(1261) . (|27p for states of arbitrary permutational symmetry. 
While p ~ N/L is the system polarization, the relation 
between the temperature T and the parameter p = r/L is 
fixed by the condition of the minimum of the free energy 
of the whole system defined by the spectrum @ and its 
degeneracy. It has the form (see [23[, [24| ) 



1 



J 



2T (1 - 2p) 



In 



l_p 

P 



(39) 



The scaling of the thermal VNE across a phase tran- 
sition, which occurs in the system with infinite range in- 
teractions at finite temperature T c ^ [U, has been 
considered in Rcf. [23|, (24[. In this case it was shown 
that the VNE of a block od size n scales as 



S(n) = 



\ log n + j log 2nepq 
nH(p) + \ log n + C{p, p) 
nH(mm(p, 1 — p)) 



where H(a) = —a log a — (1 — a) log(l 
C(p,p) does not depend on n. 



for T = 
for < T < T c 
for T>T C 

(40) 

, q = 1 — p and 
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2.5 




Figure 4: Renyi entropies for symmetric ground states of the 
Heisenberg chain versus polarization p = N/L for the cases 
a — 1,2,3,4,6,8 (curves from top to bottom, respectively) 
and for n = 20. Continuous curves refer to the analytical 
expression in (|49|) while the full dots are obtained from exact 
calculations using RDM eigenvalues in (|26l) (1271) for L = 2 * 
10 3 and same polarization. 

Another quantity of interest strictly related to the 
entanglement entropy is the mutual information, Iab, 
which measures the work necessary to erase all correla- 
tions in the bipartite system 

Iab = S(A) + S(B) - S(AB), (41) 

where S(X) is the VNE of the subsystem X. At nonzero 
temperature we find for a subsystem of size n of a system 
with size L, using (|4TJ)) : 

I AB (n, L,T) = - \og(n{L - n)) + Const (42) 

for all T <T C and I A b(T) = for T > T c . 

VI. MOMENTS OF THE REDUCED DENSITY 
MATRIX AND RENYI AND TSALLIS 
ENTROPIES 

Besides the entanglement entropy, the Renyi, R(a), 
[H and Tsallis, (T(a), [27| entropies, defined as 

b R (a, n) = — , S T (a,n) = , 

1 — a 1 — a 

(43) 

with a a positive real number, are also commonly used as 
a measure of entanglement. Notice that both expressions 
reduce the VNE in the limit a — > 1. The knowledge of 
these generalized entropies requires the computation of 
Tr(p" n - ) which, except special cases (see below) it is a very 
difficult task. For a positive integers, however, the mo- 
ments Tr(p" n s j can be computed using a quantum field 
theory (QFT) procedure which is known as the replica 
method [28| (reminiscent of the "replica trick" of disor- 
dered systems) . In this case the entanglement entropy is 



2.5 




Figure 5: Same as in Fig. U but for Tsallis entropy. Con- 
tinuous curves refer to the analytical expression in (|50[) while 
the full dots are obtained from exact calculations using RDM 
eigenvalues. 

obtained through an analytical continuation of Tr(p" n ^) 
from positive integers to real a values, using 

S { n) = -Itm^-^-Trip^f}. (44) 

In the case of 1+1 conformal field theories critical mod- 
els at zero temperature (for ground state) the displays 
universal properties, namely 

_ a-l/a 

Tr(p? n) )=C a l- S m—) (45) 

where c is the central charge of the underlying conformal 
field theory. Similarly, for the quantum XY chain with 
periodic boundary conditions at zero temperature it has 
been shown that the RDM is independent on the block 
size n and the moments can be expressed in the form (29j 

i , n°° n + P ~ 2m ^ a ) 

(^(P? n) )) i/2 = |j^ ( ( ;+ e e _ 2m£) j (46) 

where e depends on the anisotropy and transverse field 
parameters. Except these and few other cases, analyti- 
cal properties of RDM for interacting systems are largely 
unexplored. The characterization of the RDM spectrum 
given for permutational invariant systems allows to pro- 
vide another exact result for the RDM moments which 
is not accessible by QFT methods (our model is not con- 
formal invariant). 

In particular, the ground states of the ferromagnetic 
Heisenberg chain being characterized by the YT with 
r = 0, are fully symmetric states with respect to the per- 
mutations (see appendix). For these states then one can 
obtain the analytical expression of Tr(p a ) straightfor- 
wardly, using the Gaussian distribution of the symmetric 
RDM eigenvalues derived in ()35l) . The approximation 
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2~2 k ... « nj ...dx, indeed, readily provides 



VII. CONCLUSIONS 



Tr{p a ) « n 



infty 



-infty 



1 / n 2 (p-x) 2 

exp I ^—^ a | ax 



(2na 2 ) a / 2 



2<j 2 



1 i=a i 

= -= (27rpq) 2 ra~ 



(47) 



In Fig. [3] we compare the behavior of the traces of the 
RDM powers with the block size n, as obtained from Eq. 
(|47|) and from exact expressions of RDM eigenvalues. We 
see that the agreement is very good this confirming the 
correctness of our analytical derivation. 

In the case the original global state has the form of 
a maximally mixed state, i.e. is the sum of equally 
weighted projectors on symmetric states \^l,n) of the 
form p = j^t[Y,n=o \^l,n)(^l,n\, the reduced density 
matrix has one eigenvalue only Xk = ^rj, which is de- 
generate n + 1 times. In this case, then 



(48) 



Tr(pf n) ) = (n + iy 



Notice that the entanglement entropy at T = in (|40|) 
follows from (|4"7| using the expression (|44[) of the QFT 
replica method. Another quantity directly related to 
the RDM moments is the effective dimension defined as 
d e ff = T r ( p -z — j- Summarizing the above results, we have 
for this quantity that: 



To summarize, we have provided explicit analytical ex- 
pression of the reduced density matrix of a subsystem 
of arbitrary size n of a permutational invariant quan- 
tum many body system of arbitrary size L and charac- 
terized by a state of arbitrary permutational symmetry. 
We have shown, on the specific example of the spin 1/2 
Heisenberg model, that the RDM acquires a block diag- 
onal form with respect to the quantum number k fixing 
the polarization in the subsystem conservation of S z ) and 
with respect to the irreducible representations of the S n 
group. Analytical expression for the RDM elements and 
for the RDM spectrum arc derived for states of arbi- 
trary permutational symmetry and for arbitrary fillings. 
These results are provided by Eqs. ([21]), (J22J) and (|2"7|) 
presented above. Entanglement properties have been dis- 
cussed both in terms of the VNE and of the Renyi and 
Tsallis entropies. In particular, the temperature depen- 
dence and the scaling of the VNE across a finite tempera- 
ture phase transition have been considered and the RDM 
moments and the Renyi and Tsallis entropies have been 
calculated for symmetric ground states of the Heisenberg 
chain and for maximally mixed states. These results be- 
ing based only on the permutational invariance and on 
the conservation of S z (number of particles for non spin 
systems) are expected to apply also to other quantum 
many-body systems with the same symmetry properties. 



~ n for maximally mixed symmetric state 

~ y/n for pure symmetric state 
~ n x l A for critical XXZ model ground state 
n c / 4 f or a critical state with central charge c 



From the expression of Tr(p" n ^) in [47] the Renyi and 
Tsallis entropies for fully symmetric states follow as 



a i \ 1 ,„ s log a 

S R (a,n) = -log(2impq) - — -, 

2 2(1 — a) 



St(&, n) 



I (2impq) 2 — 1 

Jo. 1 — a. 



(49) 
(50) 



In Figs. [4] [5] we compare the above analytical ex- 
pressions for the Renyi and Tsallis entropies with ex- 
act calculations using the RDM eigenvalues in Eqs ([2"o]) 
(|27[) . from which we see that a very good agreement is 
found. Also notice that in the limit a = 1 both en- 
tropics reduce to the entanglement entropy (j4T)]) at T = 0: 
Sr(1, n) = = 2 1°S ^Trepqn. 

In general, for arbitrary permutational symmetries and 
for finite temperatures, one must recourse to direct cal- 
culations using the general expression ([2"o]l for the RDM 
eigenvalues, since it is not easy in these cases to give 
simple analytical expressions of a. The study of the ana- 
lytical properties of the RDM moments represents an in- 
teresting problem which deserves further investigations. 
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Appendix A: RDM elements for symmetric states 

In this appendix we provide a proof of Eq. (|21[) which 
is valid for fully symmetric states (case r = 0) such 
as, for example, the ground state of the ferromagnetic 
Heisenberg chain. For r = 0, the corresponding YT is 
nondegenerate and the state of the full system is pure: 
P = \^l,n)(^l,n\ with \^l,n) the symmetric state 



I* 



L.N 



-1/2 



(Al) 



where the sum is over all possible permutations. Since 
all L sites are equivalent due to permutational invariance, 
any choice of n sites i\, 12, i n within L sites gives the 
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same RDM, which we denote by p^ 



L,N,0 
n) 



Trt-nP- It 



has been shown in [10( that p^^'° takes form 



P 



L,N,0 
(n) 



n f L—n\ 



In the natural basis the matrix elements of RDM are 
given by the above discussed values of observables. Using 
(|Alj) . one explicitly computes all RDM elements as 



L,N,oY 3 - m 



( L,N,0\ ( L.N,o\ 



£ i +i + ■ ■ + ■ m \N—% 
...+m' / L\ 



. ■ 

(A2) 

with w = i + j + ... + m (the sets ij...m and i'j'...m' 
are binary representation of numbers P — 1, Q — 1). In 
this case we have that all elements of a block Bk are 
equal (this is not true for r > 0). We also see that the 
the elements in Eq. (|A2[) are the same as those obtained 
from Eq. (f2Tj) for r = 0. Note that in the thermodynamic 



limit r\ = 1, and 
with Eq. 



IL—n 
\N~k 



)/( 



! (l—p) in agreement 



Appendix B: RDM elements in the thermodynamic 

limit 



To calculate the RDM, we shall use the representation 
([7p for the density matrix of the whole system tr, rewrit- 
ten in the form 



P(n) 




Trr 



Ttl- 



(Bl) 



Note that the L\ permutations can be done in three steps: 
first, choose at random n sites i\ ^=12 7^ ■■■ 7^ in among 
the L sites. There are (^) such choices. Then, permute 
the chosen n sites, the total number of such permutations 
being n\. Finally, permute the remaining L — n sites, the 
total number of such permutations being (L — n)!. The 
latter step (c) under the trace operation is irrelevant be- 
cause these degrees of freedom will be traced out. The 
operation permuting n sites commutes with the trace op- 
eration since Tr^- n does not touch the respective subset 
of n sites. Consequently, (|B1[) can be rewritten as 



(B2) 

We recall here that a filled YT of type {L — r, r}^ con- 
tains a mixed symmetry part with 2r sites and r spin 



up , and a fully symmetric part with L — 2r sites and 
N — r spin up (in the following we adopt an equivalent 
terminology which refers to spins up as particles and to a 
spins down as holes). This implies that the correspond- 
ing wave function \^L,N,r) factorizes into symmetric and 
antisymmetric parts as 

* = |012)®|^34)®---®|02r-l,2r)®|*L-2r,AT-r)2r+l,2r+2,..L 

(B3) 

with the antisymmetric part consisting of the first r fac- 
tors of the type 



|0 13 > = ^=(|10>ia - |01) 12 ) 



(B4) 



and with the symmetric part, |\& L-2r,N—r)i given by 
(|A1[) . A general property of factorizcd states implies that 
if the global wave function is factorized, |$) = \ip)i\4>)n 
and out of n sites of the subsystem, m sites belong to 
subset /, and the remaining 712 = n — n\ sites belong to 
the subset iT, then the reduced density matrix factorizes 
as well: 



(B5) 



To do the averaging, we note that among total number 
of choices (^) there are (a) ( L ~ 2r ) possibilities to choose 
n sites inside the symmetric part of the tableau, contain- 
ing N — r particles, (b) 2r(^ i ~ 2 ] r ) possibilities to choose 
n — 1 sites inside the symmetric part of the tableau and 
one site in the antisymmetric part (c) ( 2 2 r ) 2r ) possi- 
bilities to choose n — 2 sites inside the symmetric part 
of the tableau and two sites in the antisymmetric part 
and so on. The contributions given by (a) and (b) to the 
right hand side of (|B2|) for ' r are given, according to 
m, by 



F 



F,M,0 
P(n) ' 



2r 



F,M,0 



(B6) 



with F = L - 2r, M = N - r and with pi = 1/2 the 
density matrix corresponding to a single site in the an- 
tisymmetric part of the tableau. Brackets (.) = ^ ^ p 
denote the average with respect to permutations of n ele- 
ments. The contribution due to (c) to (|B2[) splits into two 
parts since the ( 2 2 r ) possibilities to choose two sites in the 
antisymmetric part of the tableau consist of 4 ( 2 ) choices 
with two sites into different columns and the remaining 
r choices with both sites belonging to a same column. 
For the former choice, the corresponding density matrix 
pi , while for the latter case is given by 



F.M.O 

18 P(n-2) » P\ 

F.M.O „ 1 
P(n-2) ® Pasymrn, With 



l_L (1 „-oi))(-L ( io-oi)| = I 



/0 ON 
1-10 
0-110 

Vo 00/ 

(B7) 
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Proceeding in the same manner for arbitrary partitions 
of Z sites in the antisymmetric part of the tableau and 
n — Z sites in the symmetric part, we get 



L.N.r ( L 



§Q Pasymm J 2 




F \ F.Mfi \^ i 
_ Z )P(n-Z) 2^ 




(B8) 



formula is only true for k = Z, otherwise elements con- 
stituting Gz(W) are not all equal). Restricting to the 

case k = Z and denoting W m = p^f^ ® (p± 
have 



, we 



in—m\ 

r(w m ) = r m = ±^1. 



(B12) 



\2Z1 



It is worth to note that the element gz C Gz is simply 
given by 



From this the general scheme for the decomposition of 
the general RDM becomes evident. In the above formula, 
Q 

the products Y[ with Q < i are discarded. The matrix 

i 

elements pf^\'° ar c given by (|A2I) . 

For simplicity of presentation, we prove Eq. (|22[) for 
the case Z = k and then outline the proof for arbitrary 
Z. 

In the thermodynamic limit one can neglect the differ- 
ence between factors like 4(£) and ( 2 2 r ) in Eq. (|B8|) . The 
latter then can be then rewritten in a simpler form as 



L 



L.N.r 
P{n) 



F 



F,M,0 , I 2r W r \ F ,M,0 

n/'W f 1 U-l Vi)' 



F 



(B9) 



2r 

2 )\n-2 



F \ F.M.O „ 

P(„_ 2) ®pi O pi 



Note that one can omit all terms in (|B8[) containing 
Pasymm since the respective coefficients correspond to 
probabilities of finding two adjacent sites in the asym- 
metric part of the YT (proportional to r), which vanish 
in the thermodynamic limit, respect to the total num- 
ber of choices which is of order of r 2 . A sub-block Gz 
of a block k consists of all elements of the matrix pi n ) 
having Z pairs of e^ej 1 in its tensor representation, like 
e.g. 



-1Z 



Z + 21 + 12 + ■■■ + i n -2Z = k. The total number of elements 
gz C Gz in pf'^' r is equal to the number of distributions 

of Z objects el, Z objects ej, and (k — Z) objects e\ on 
n places, given by 



deg G z 



Z\Z\{k - Z)\(n - k - Z)\ 



(BIO) 



(this is another way of writing (fT4"]) V Each term W in 
the sum (|B9|) after averaging will acquire the factor 



T(W) 



deg G Z (W) 
degGz 



(Bll) 



where deg Gz{W) is a total number of gz elements in 
the term W, provided all of them are equal. For in- 
stance, degG z (p™'°) = degGz, deg G z (p^-m) ® 

(PiY" 1 ) = ( 2 #)r 2 T)EZ=o(Z)(rf-Z)( last 



>-r.Qrf- + r 1 (t)( I1 ! 1 ) l£ ^- 



2r\( F \9o 



(n-2,k) 



2 \n-2 



with q=l-p and g { ^ k) - - 



(B13) 



1-2/i 



(l-2^) 



is the element of a p^ 1 ' corresponding to a block with 
k particles (the factors T m are due to the averaging while 

the factors come from (pi ) )■ Restricting to the 
case k = Z, and taking into account 

L-m) „ n! (1 - 2 p)" ^ _ (2/x) 
(£) ~(n-m)! F" 

so that 



L ) \m 2 m I my 



(B14) 



we finally obtain, using (|B13[) . that 



%-2Z 

- E /*" 

ro=0 



n-2Z 



= (p-/i) n - z (g-/i) z £ 



p m /n - 2Z 



(p — p) m \ m 



(p-p) n - Z (q-p) Z ! 



p n-Z q Z f (P~P)(<1-P) 

V PQ 
p n - z q z V z = V z g , 



P - P 
z 



n-2Z 



(B15) 



with go the diagonal element in the same block k = Z . 



In the last calculation we used the relation 



( n ~m. Z } • This proves formula (f2"2"j) for the particular case 
k = Z and arbitrary n. 

For arbitrary k, Z, one proceeds in similar manner as 
for the case k = Z case. Since the respective calculations 
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are tedious and not particularly illuminating, we omit 
them and give only the final result: 

g z = » m E (p-») n - m - k+i 

m=0 i=max(0,Z+k — n+m) 

( q -,) k -f- Z )( n - k -. Z ), (B16) 
\ i J \ m — i J 

which, after some algebraic manipulation, can be rewrit- 
ten in the form 

«%{^t, Z )=^'-^ (BIT) 



This concludes the proof of Eq. pip in the thermody- 
namic limit L — > oo. 
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